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ABSTRACT 



We present a versatile numerical algorithm for computing resonances of open dielectric cavities. The emphasis is 
on the generality of the system's configuration, i.e. the geometry of the (main) cavity (and possible inclusions) 
and the internal and external dielectric media (homogeneous and inhomogeneous). The method is based on a 
scattering formalism to obtain the position and width of the (quasi)-eigenmodes. The core of the method lies 
in the scattering S-matrix and its associated delay Q-matrix which contain all the relevant information of the 
corresponding scattering experiment. For instance, the electromagnetic near- and far-fields are readily extracted. 
The flexibility of the propagation method is displayed for a selected system. 



1. INTRODUCTION 

In the research areas of microcavity physics [ 1 1, lasing [2] and sensing [3 1, increasing flexibility in the geometries 
of the cavities and in the constituting/surrounding media has been paramount for better performance regarding for 
instance, directional laser emission, lower lasing threshold and higher sensing resolution. Confronted with this 
demand, we have set out to develop a versatile theoretical approach capable to meet the challenges of ever more 
complex experimental arrangements. Although finite element methods offered themselves as a possible solution, 
we thought that they somewhat mask the underlying physics. A scattering approach appeared to us as the method 
of choice complementary to existing boundary element method (BEM) [4| and boundary integral method (BIM) 
to name a few. We have extended and upgraded a method |6 7] recently introduced that has the potential to 
treat a wide variety of geometrical and refractive index deformations. This presentation is but a glimpse at the 
capabilities of the method. Details will appear elsewhere. 



2. TRANSFER, SCATTERING AND DELAY MATRICES 

Our tale is that of a set of matrices: transfer matrices to propagate the solution from the interior of a system to 
the outside world, a scattering S-matrix to collect the information and a derived delay Q-matrix to provide the 
eigen-modes of the system and their quality factors (or delays). There seems to have been renewed interest of late 
on the properties of the Q matrix in different contexts, and we refer the reader to some recent publications for 
further details EH3E3HI1. 

Here is a brief summary of the important elements of the method. For the sake of simplicity, let us consider the 
solution \|/(r) of an electromagnetic wave in a two-dimensional (2D) dielectric medium. We will concentrate on 
TM polarisation. The equation to solve is the Helmholtz equation 



[V 2 + n 2 {r)k 2 ] y(r) = 



(1) 



which, in polar coordinates, takes the form 




+ 



r^ + -^+n 2 (r^)k 2 r 2 y(rrf)=0. 



(2) 



We then proceed to divide space into N +1 interfaces (J = 0,1,... ,N) located at distances Rj = ro + 2je (Fig. 
1(a) l. To follow the procedure, let us focus our attention on a thin annular region of thickness 2e and central 



(a) Concentric domains 



(b) Expansion coefficients 



Figure 1: (a) Schematic representation of the separation of a generic cavity in annular domains. The refractive 
indices, n; n and n„, are constant within the innermost disk and outside of the outermost domain respectively. \(b) 
The expansion coefficients a' and b J at the successive concentric domains are used to propagate the solution to 
the exterior region where the S matrix (R = SAj is obtained. 



radius rj, j = 1,2. ..N, covering part of the generic cavity. Over this region, we assume the refractive inde^to be 
dependent upon (j) only and n 2 (r,§)r 2 is then evaluated at rj for all (]). Turning to (|2j, we obtain an approximate 
local expansion of the differential equation over \r— rj\ < £, 
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V(r,4>)=0. 



(3) 



This differential equation is separable for the local wavefunction \\t J \r, (j)) . Under the Ansatz \|/ ; (r, 0) = %] (r)<t>i ( 
we obtain two exact differential equations 
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ZJiPj) = o 



<t> J 



= 



(4) 
(5) 



where pj — r/rj and <tj is a separation constant. Imposing periodic boundary conditions on Eq. (|5jl and expanding 
4> ; ((j)) in a Fourier series, we end up with an eigenvalue problem of an hermitian symmetric matrix. The re- 
trieved eigenvalues specify the set of constants {^}, and the associated eigenvectors define a basis of normalized 
eigenfunctions {^((j))}. We may then solve the Cauchy-Euler differential equation Q exactly, 



(6) 



where ajj and bjj are constants to be obtained from the boundary conditions. The sign of ^ renders ^/(p,-) 
evanescent (- sign) or propagating (+ sign) (Fig. ( 1(b) I. 

This procedure is propagated to an adjacent thin annular domain of thickness 2e centered on r, + i = r, +2e. 
Boundary conditions for TM polarisation assume continuity of and and their normal (radial) derivative at 
their common boundary p /+ = 1+z/rj and p 7 + 1 - = 1 — e/ rj + 1 . 



Finally, having reached the outermost region (r > R„ 
incoming and outgoing H^ + ' Hankel functions 



rw + e), the solution is matched to a superposition of 



y(r,^)=£ A„,H,n } {n kr)+B m H£ } {n kr) 



(7) 



1 The refractive index is real in the present description. However, the method can easily accommodate a complex refractive index to account 
for absorption and loss processes. 



from which the scattering S-matrix is defined 

B = SA. (8) 

One can then show, and this is a crucial part of the approach, that the characteristic expansion vectors A are nothing 
less than the eigenvectors A p of the delay Q-matrix lfTTl[T2 l given by 

such that 

QA p = cx p A p (10) 

where cz p (in length units) are the eigen-delays of a characteristic mode p proportional to its quality factor or 
inversely proportional to its (quasi)-resonance width. 

Clearly, this short description does not do honour to the daunting task of implementing the propagation steps, 
transferring information between adjacent annular domains, nor does it do justice to the theoretical subtleties 
necessary to interpret the results obtained. However, it should suffice to give the jest of the method. An example 
will give a glimpse of its capabilities.. 



3. RESULTS: A SELECTED EXAMPLE 



Figure (2(a) i present some of the geometries that we have used for the calibration of the method: on the top 
row, we see a number of homogeneous cavities of increasing distortions with respect to a perfect disk, while the 
discontinuities in the dielectric media increase from top to bottom. We present in Figs ( |2(b)| - j2(c)| i results of a 
calculation of a more complex arrangement with both a geometrical distortion and a dielectric inclusion (in form 



of a C, for complex no less). The sweep over a wavenumber window (Fig. 2(b) i reveals a number of resonances of 



different quality factors: the ordinate axis displays the eigenvalues of the Q matrix and are directly proportional to 



the width of the corresponding resonances. In Fig. (2(c) i, we show the near- and far-field distribution for 3 selected 
resonances appearing in the spectrum. They have been chosen for their respective differences. Resonance (A) has 
the near-field (NF) structure of a whispering gallery type, almost uniform, whereas its far-field (FF) distribution 
shows pronounced peaks at certain privileged directions. Resonance (B) has a NF amplitude completely located 
within the inclusion and a FF with a dominant distribution in a direction perpendicular to the opening of the 
inclusion. Resonance (C), a low quality mode, has for its part strongly non-uniform NF and FF distributions. 

Typically, our method handles well complicated combinations of various geometries and medium discontinuities. 
However, it is fair to say that its careful implementation has only been thoroughly tested for two-dimensional 
systems. Among other things, the three-dimensional cases are waiting in the wings ... 



4. CONCLUSION 



We have presented a scattering/propagation method applicable to cavities of arbitrary shape and arbitrary inhomo- 
geneities of the medium (continuous and/or discontinuous). The calculational approach has now reached maturity 
for 2D systems, and is currently being extended to photonic molecules (combination of a few cavities), and more 
generally, to photonic complexes (periodic or aperiodic arrangements of cavities). Details of the implementation, 
and further results will be published elsewhere. 
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Figure 2: (a) Typical arrangements of geometrical distortions and/or inhomogeneous dielectric media, the com- 
bined complexity of the systems increase along the diagonal. \(b)\ Representative eigen-delay spectrum for the C 
inclusion, (c) Intensity distribution of the near-field (NF) and the far-field (FF) in false colour (the darker the 



colour, the higher the intensity). The outer-circle serves as a support to display the FF. 
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